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Abstract 

It is known that if the dimension is a perfect square the Clifford group can be 
represented by monomial matrices. Another way of expressing this result is to 
say that when the dimension is a perfect square the standard representation of 
the Clifford group has a system of imprimitivity consisting of one dimensional 
subspaces. We generalize this result to the case of an arbitrary dimension. Let 
k be the square-free part of the dimension. Then we show that the standard 
representation of the Clifford group has a system of imprimitivity consisting 
of fc-dimensional subspaces. To illustrate the use of this result we apply it to 
the calculation of SIC-POVMs (symmetric informationally complete positive 
operator valued measures), constructing exact solutions in dimensions 8 (hand- 
calculation) as well as 12 and 28 (machine-calculation). 
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1. Introduction 

The Clifford groupQ plays a major role in many areas of quantum information. 
When expressed in terms of the standard basis many of the Clifford unitaries are 
Hadamard matrices, with the property that every matrix element has the same 
absolute value = 1 / \/N (where N is the dimension) . It is therefore somewhat 
remarkable that, in the special case when the dimension is a perfect square, a 
simple change of basis will cause every Clifford unitary to become what could be 
regarded as the opposite of a Hadamard matrix, namely a monomial matrix, with 
only one non-zero element in each row and each column [5] . In the following we will 
generalize this result to arbitrary dimensions. Given any dimension N there exist 
unique integers k, n such that k is square-free (i.e. k is not divisible by the square 
of any prime number) and N = kn 2 . We will show that the Clifford group admits 
a representation entirely in terms of matrices with the property that, when written 
out in block form, each row and each column contains exactly one non-zero k x k 
block. We refer to such matrices as fc-nomial. It will be seen that the result proved 
in Ref. [5] is simply the special case of this for k = 1. 

This result can be re-phrased. Let H be JV-dimensional Hilbert space. Then our 
result is that % can be written as a direct sum 

n-l 

n = H r , s (1) 

r,s=0 

where the subspaces T-L r ,s are all /c-dimensional and where, for every Clifford unitary 
U and every pair r, s 

UHr^s = Tir'.s' (2) 

for some r',s'. Such a set of subspaces is called a system of imprimitivity [H[7]. 
Systems of imprimitivity play an important role in the theory of group represen- 
tations as they can be used to construct an induced representation for the whole 
group starting from a representation of a subgroup. No such application is in ques- 
tion here since we already have a representation for the group. There are, however, 
other applications, as we will see. 

Many interesting physical properties are revealed by a deeper understanding of 
the relevant group structures and properties. The Clifford group in its various forms 
is a good example of this, having applications to fault tolerant quantum comput- 
ing [51(5] , measurement based quantum computing [TOlfTTj , describing the boundary 
between classically simulable and universal quantum computing [IHEI], quantum 
non-locality [2], mutually unbiased bases (or MUBs) [2lHl [TSH2u] . and symmetric 
informationally complete measurements (commonly referred to as SIC-POVMs or, 
as here, simply SICs) [4ll5l ll8l[l9ll2TH42] . Although the result proved in this pa- 
per is potentially relevant to all these areas we will illustrate its use in connection 
with the SIC existence problem, as that is the route by which we were led to it. 
SICs are important practically, with applications to quantum tomography [43U46j . 
quantum cryptography pTHSTj . quantum communication [52-56 , Kochen-Specker 
arguments [57], high precision radar [58H60] and speech recognition [61] (the last 
two applications being classical). They are also important from a mathematical 

■'"We should remark that there are several versions of the Clifford group [TH4], In this paper 
the focus will be on the version defined to be the normalizer of a single copy of the non-Galoisian 
version of the Weyl-Heisenberg group. However, we believe our results are likely to be relevant to 
the other versions. 
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point of view, as giving insight into the geometrical structure of quantum state 
space [62) , and from a foundational point of view, as playing a central role in the 
qbist approach to the interpretation of quantum mechanics [63H65] . SICs have been 
realized experimentally [50,66]. They have been calculated numerically for every 
dimension < 67 and exact solutions have been constructed for dimensions 2-16, 19, 
24, 35 and 48 (see Ref. [36] for a comprehensive listing of solutions known in 2010 
and Ref. [5] for the exact solution in dimension 16). There are therefore grounds for 
conjecturing that SICs actually exist in every finite dimension. However, in spite 
of strenuous efforts by many investigators over a period of more than 10 years the 
question is still undecided. The connection with the Clifford group comes from the 
fact that the overwhelming majority of known SICs are covariant with respect to the 
Weyl-Heisenberg group (or generalized Pauli group). Moreover, Weyl-Heisenberg 
covariant SICs have been found in every dimension for which SICs have been found 
at all. It is therefore tempting to conjecture that, not just SICs, but specifically 
Weyl-Heisenberg covariant SICs exist in every finite dimension. 

The investigation reported here was motivated by two striking facts concern- 
ing Weyl-Heisenberg covariant SICs. Firstly, the defining equations are massively 
over-determined. Specifically, one has d 2 equations for only 2d — 1 real param- 
eters. So if one did not know better one might think it unlikely that solutions 
existed at all, in any dimension > 2. The fact that solutions actually do exist (at 
least numerically) in every dimension up to 67, suggests the presence of some kind 
of concealed symmetry or other algebraic feature of the equations which forces a 
solution notwithstanding the over-determination; and it further suggests that iden- 
tifying that feature might be the key to the existence problem. The second striking 
fact is that the known exact solutions all turn out to be expressible in terms of 
radicals. Again, this could not be anticipated in advance of actually calculating 
the solutions, and seeing that they are so expressible. The defining equations of a 
Weyl-Heisenberg covariant SIC fiducial vector are quartic polynomials in the real 
and imaginary parts of the components. The standard way to solve such a system 
of equations (except in the simplest cases) is to construct a Grobner basis [67j . 
thereby reducing the problem to that of solving polynomial equations in a single 
variable. It turns out 24,26,31,36 that the single variable equations are all of 
degree 3> 4 and, as is well-known, such equations are only solvable in terms of 
radicals in the exceptional case when the Galois group is solvable [69]. In fact, 

one findU [M2H2] that 

in all the known cases the Galois group has a particularly 
simple form even amongst the class of solvable groups, having a subnormal series of 
the form {e} <3 H <\ G where e is the identity, H and G/H are Abelian, and G/H 
is order 2 [G being the Galois group of the smallest normal extension of Q con- 
taining all the components of the fiducial together with \/— 1). It seems probable 
that this feature of the Galois group is giving us an important clue, and we would 
naturally like to investigate it further. The trouble is that, although the Galois 
group has a very simple structure, the solutions themselves are not (in many cases 
the print out extends to many pages), which tends to obstruct progress. We were 
therefore led to wonder if a change of basis would have the effect of simplifying the 
solutions. In Ref. [5] it was shown that the monomial representation does indeed 



2 As originally noticed by Jon Yard | 68| . although Ref. 1301 gives correct expressions for the 
exact fiducials in dimension 14, there is an error in the calculation of the Galois group. The 
correct Galois group for this dimension can be found in Ref. |42| . 
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produce considerable simplification for N = 4, 9, 1 6. In the following we show that 
the fc-nomial representation is similarly efficacious. We perform a hand calculation 
of the Weyl-Heisenberg fiducials for N = 8 and we find that the expressions which 
result are indeed much simpler than the standard basis expressions obtained by 
machine calculation |36j . The fact that a hand calculation is even possible may be 
taken as further evidence that the fc-nomial basis is better adapted to the problem. 
We also describe a machine calculation of a fiducial in dimension 28. Although 
exact fiducials have previously been calculated in dimensions 24, 35 and 48, these 
relied on the existence of additional symmetries of order 6, 12 and 24, respectively. 
So far, N = 28 is the largest dimension in which it has been possible to find an 
exact solution only with the help of a symmetry of order 6, while N = 16 in Ref. [5] 
is the largest dimension for which the computation relied only on a symmetry of 
order three, as conjectured by Zauner. This may be taken to confirm the hypothesis 
that the fc-nomial basis is indeed better adapted to the SIC problem. As further 
confirmation of this point we revisited the calculation of an exact fiducial in di- 
mension 12 first reported in Ref. [31] , We found that the use of a fc-nomial basis 
reduced the computation time by more than 3 orders of magnitude. We also found 
that the expressions for the fiducial are much more compact than the ones given in 
Refs. [3TJISS] 

Finally, let us observe that the advantage of a hand calculation, such as the 
solution for N — 8 presented here, or the solution for N = 9 presented in Ref. [5], 
is that it gives us a degree of insight into the algebraic intricacies of the problem 
which a machine calculation cannot provide. We suggested above that the fact that 
solutions exist notwithstanding the massive degree of over-determination, and the 
striking simplicity of the Galois group, both hint at the presence of some underlying 
symmetry, or other algebraic feature of the equations which has so far eluded us. It 
may be that a study of hand-constructed solutions will lead us to the secret which 
will enable us finally to crack the problem. 

The plan of the paper is as follows. In Section [2] we describe the fc-nomial 
representation of the Clifford group. In Section [3] we briefly review a few basic 
facts about SICs. In Section [4] give a hand calculation of the Weyl-Heisenberg SIC 
fiducials in dimension 8 . In Section [5] we obtain, by machine calculation, a very 
compact solution for dimension 12. In Section|6]we obtain, by machine calculation, 
exact fiducials in dimension 28. 



To fix notations we begin by describing the standard basis representation of the 
Clifford group. For more details see (for example) Ref. [25]. Let |0), . . . \N — 1) be 
the standard basis in dimension N. Define the operators X and Z by 



where addition of ket-labels is mod N and u) — e 2 «* . The Weyl-Heisenberg dis- 
placement operators are then defined by 



2. The fc-NOMiAL Representations 



X\u) = \u+l) 
Z\u) =gj u \u) 



(3) 
(4) 



D p = T PlP2 X Pl Z P2 



(5) 
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where p = ( p 2 ), t — — e « and p\, P2 run from to N — 1, where 

N={ N N ° M (6) 

[2iV N even 

With this definition 

£> p £> q = r<™>D p+q (7) 
where the symplectic form (•, ■) is defined by 

(P,q) =P29i -i>i92 (8) 
The symplectic group SL(2, Z^y) consists of all matrices 

such that a, /3, 7, 5 G Zjy and DetF = 1 (mod TV). To each such matrix there 
corresponds a unitary L/p, unique up to a phase, such that 

U F D p Ul = D Fp (10) 

If /3 is relatively prime to N we have the explicit formula 

18 N-l 

u F = -j= 2^ T y >m{v\ (11) 

where /3 _1 is the multiplicative inverse of (3 (mod iV) and e %e is an arbitrary phase. 
Observe that it follows from this formula that Uf is a Hadamard matrix in the 
standard basis whenever (3 is relatively prime to N. If (3 is not relatively prime to 
N we use the fact [25] that F has the decomposition 



<i2) 

where /?i , /?2 both are relatively prime to N so that {/ ^ , Uf 2 can be calculated 
using Eq. (JTTJ) . £/f is then given by 

U f = U Fi Uf 2 (13) 

up to an arbitrary phase. The Clifford group then consists of all products D p IIf- 
The extended Clifford group is also important. Enlarge SL(2,Z 7 y) to the group 
ESL(2,Zjy) by including the anti-symplectic matrices 



v 7 ^ j 

with Deti* 1 = — 1 (mod N). Each such matrix has the unique decomposition 

F = FJ (15) 

where F is symplectic and 



-1 / 

To J we associate the anti-unitary Uj which acts by complex conjugation in the 
standard basis: 

(JV-l \ N-l 

u=0 / u=0 



for all \ip) = ^2^=0 i'v-Wl- We then associate to F the anti-unitary 



U F = U p Uj (18) 

'■<•■<-!'< 1'" n pUp, W1U1 p C ^ 



The extended Clifford group consists of all products of the form D p U F , with p G Z 2 ^ 



and F G ESL(2,Zat). 

These preliminaries completed we now turn to the construction of the fc-nomial 
representations. Let TV = kn 2 (one would usually choose k to be square- free, so 
as to make the non-zero blocks in the representation matrices as small as possible; 
however this is not necessary). Then it can be seen from Eq. (0) that X kn and Z kn 
commute, and are therefore simultaneously diagonalizable. It is easily verified that 
a joint eigenbasis is 



n-l 

\r, a, j) = -= V A- rf | (s + jn) + tkn) (19) 



where A = e »' , where r, s G Z n label the eigenspaces according to 

X fe "|r, S ,j)=A r |r, S ,i) ^ fc "|r, S ,j) = X s \r,s,j) (20) 

and where j G Z& labels the basis vectors within each eigenspace. The action of X, 
Z in this basis is 

!\r,s + l,j) s^n-1 
\r,s + 1J+1) s = n-l,j^k-l (21) 

\ r \r,s + l,j + 1) s = n-l,j = k-l 

Z\r,sJ)=u s+n i\r-l,s,j) (22) 

from which we see that the displacement operators are fc-nomial in this basis (in 
fact monomial) just as they are in the standard basis. From Eq. (|2"2")l it follows that 
not only X kn and Z kn are diagonal in this basis, but also Z n . 

The advantage of this basis is that in it, not only the displacement operators, 
but also the symplectic unitaries are represented by fc-nomial matrices. In fact, let 

F=t f) (23) 



7 5 

be an arbitrary matrix G SL(2,Z^). Then 



U F X^U F = I i-lV^Z-^ k odd, n even 

1 X kn8 z -kn 7 otherwise 



UlZ^U F = [ tT*:!T Zkna k ^;Z^ (25) 



x -knp z kn a otherwise 



where we used the fact that 

_fe 




N odd 
TV even 



(26) 



to calculate the signs. Consequently 

X kn U F \r,s,j) = \ r 'u F \r,s,j) (27) 
Z kn U F \r, S ,j) = \ s 'u F \r, S ,j) (28) 







where 

r = Sr — 7s + mjS (29) 
s' = -f3r + as + ma/3 (30) 

with 

f§ k odd, n even, 

m = < z (31) 
I otherwise. 

So Uf takes the (r, s) eigenspace to the (r',s') eigenspace: 



fc-i 

U F \r, s,j) = Yl ( M ^rs) jf \r', s',f) (32) 

for some family of k x fc matrices M F ,rs- It is thus fc-nomial, as claimed. Note that 
this also shows that the eigenspaces labeled by (r, s) form a system of imprimitivity 
(see Eq. |2])). 

To calculate the matrices Mp rs suppose, first, that /3 is relatively prime to N. 
Then using Eq. (1111) one finds, after a certain amount of algebra, that, up to an 
overall phase, 

{M F , rs )-, = —r^W+y^-^'+y^+^+^+^f) (33) 
n \fk 

where s' is given in terms of r, s by Eq. (|30[) (note that in this formula it is essential 
that s' be reduced mod n, so that it is an integer in the range < s' < n). The case 
when /3 is not relatively prime to N can then be handled using the decomposition 
of Eq. p]). 

It is perhaps interesting to note that when (3 is relatively prime to N the non-zero 
blocks in the fc-nomial representation are Hadamard (in contrast to the standard 
representation where the whole matrix is Hadamard). 

Finally, we can calculate the anti-symplectic anti-unitaries using Eq. (fT8|) to- 
gether with the fact 

Uj\r )S) j) = \-r,s,j) (34) 
3. SICs 

The A:-nomial representations described in the last section are a very general con- 
struction which is potentially relevant to any area where the Clifford and extended 
Clifford groups play a role. In the remainder of this paper we illustrate their use 
by showing how they can be used to simplify the calculation of SICs. We begin, in 
this section, by reviewing a few basic facts (for more details see Refs. [22 | [23 j [25 l l36] 
)• 

A SIC in dimension N is a family of TV 2 operators 

E r = ^\A){A\ (35) 

such that the vectors \ip r ) satisfy 

KlMlMl = I 1 r r " S ' . (36) 

ITf+T otherwise. 
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The fact that the matrix Tr (E r E s ) is non-singular means that the E r are a basis 
for operator space and 

N 2 

Y,E r =L (37) 
r=l 

So the E r form an informationally complete POVM (positive operator valued mea- 
sure). 

The vast majority of SICs which have been constructed to date are covariant 
under the action of the Weyl-Heisenberg group. That is, they have the property 

D^Dl = £ p+q (38) 

for all p, q (where, instead of labelling the POVM elements by the single index r, 
we label them by p £ For a SIC of this kind we can generate the entire SIC by 

acting on Eq with the displacement operators D p . The vector l^o) corresponding 
to Eo is called the fiducial vector. We usually omit the subscript and simply denote 
the fiducial vector \ip). The problem of constructing a Weyl-Heisenberg covariant 
SIC thus reduces to the problem of finding a single vector \ip) such that 

i , ( . f 1 p = mod N , 

\(1>D P 1> = i , n AN 39 

It is a striking, so far unexplained fact that every known Weyl-Heisenberg SIC 
fiducial vector is an eigenvector of a Clifford unitary D^Uf for which TrF = — 1 
mod N and F ^ I (note that the requirement F ^ I is automatic if TV ^ 3). 
It can be shown [22,25,36 that such unitaries are always order 3 up to a phase. 
Conversely, it appears that every such unitary has a Weyl-Heisenberg SIC fiducial 
as one of its eigenvectors. In the following sections we exploit this fact by looking 
for fiducials which are eigenvectors of Uf, , where 

* - (! :!) w) 

is the Zauner matrix [22]. We refer to Uf, as the Zauner unitary. 

4. Dimension 8 fiducial 

We now perform a hand-calculation of the fiducials in dimension 8 which are 
covariant under the version of the Weyl-Heisenberg group considered in this paper 
(for fiducials covariant with respect to a different version of the group see Refs. [EJ 
[2Tj). Of course exact expressions for the standard basis components have already 
been found by Scott and Grassl |36j (by means of a machine-calculation of a Grobner 
basis). So if we only wanted to make the point that the expressions become much 
simpler in the fc-nomial basis it would be enough just to transform the expressions 
in Ref. [36]. However, as we stated in the Introduction, we feel that a hand- 
calculation gives a greater degree of insight into the problem than one gets from a 
machine-calculation. In particular we entertain the hope that a close inspection of 
hand-calculated solutions may enable us to spot the "secret ingredient" responsible 
for the fact that SICs exist (at least numerically for dimensions < 67) in spite of 
the massive degree of overdetermination, and for the striking features of the Galois 
group — thereby, perhaps, leading to a solution of the existence problem. 

It may be worth remarking that we found it surprising that we were able to solve 
the equations at all (by hand, that is) . One of us (DMA) well recalls that he could 
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only obtain the dimension 7 fiducials in Ref. [25] by a process of trial and error, 
and that in order to be sure that he had found all the solutions he had to appeal 
to the numerical work of Renes et al. |23| . 

We look for SIC fiducial vectors which are eigenvectors of the Zauner unitary, 
Up, ■ Note that 8 is one of the special dimensions [36] for which fiducials can be 
found in all three eigenspaces of Uf z ■ The first step in the calculation is to choose 
a basis which diagonalizes Uf z ■ Since there are only 3 eigenspaces there are many 
ways to do this, which are far from being equivalent from the point of view of 
calculational complexity. One natural way to proceed, which causes the equations 
to be solved to take a particularly simple form, is to consider the normalizer Mz 
of the cyclic subgroup generated by Fz — i.e. the set of matrices G <G ESL(2,Zjy) 
such that 

GF Z = F r z G (41) 
for r = 1 or 2. The significance of Mz is that the corresponding unitaries and 
anti-unitaries permute the eigenvectors of Uf z - In particular they permute the 
fiducials which are eigenvectors of Uf z - Having constructed Mz one then picks 
out a maximal subset containing Fz such that the corresponding unitaries/anti- 
unitarics commute. The joint eigenvectors of the operators in this subgroup give 
us our desired basis. As we will see the fact that the matrices are all fc-nomial 
considerably facilitates the calculation. 

The group Mz is generated by the matrices 

OA , ( 1 5\ „ f-l , (42) 



K = 



1 



A = 



P = 



-1 



v -5 6 

Here K is order 2 anti-symplectic, A is order 24 anti-symplectic and P is order 2 
symplectic. We have 

F z = A 8 (43) 

and 

PA = AP PK = KP AK = KA~ l P (44) 

Define symplectic matrices 



K = 



-1 

1 



A = 



1 -5 

-5 -6 



We then have 



K = KJ 



A = AJ 



and, consequently, 



U K 



U a 



U A Uj 



(45) 

(46) 
(47) 



where J is the matrix defined in Eq. ()16[) . 

We now calculate the matrices corresponding to these operators. For dimension 
8 we have k — n — 2. If we order the basis vectors \r, s,j) lexicographically 

|0,0,0},|0,0,1),|0,1,0),...,|1,1,1) (48) 

the matrices become (with a given choice for the arbitrary phase factors) 



(K 








°\ 








Ki 








K 2 








\o 








K 3 J 



( Aq 








°^ 








Ai 














A 2 


\o 


A, 








(49) 



9 



and 



Up 



where a x , a z are the Pauli matrices and 



1 








^ 
























V° 








-CF X ) 



K 
K 3 



1 



-(' 

i A i 



i - 

-i - 

i i 



V2 v 
V2\i 



A 2 = ^ 



A, 



V2 



V2 



We fix the arbitrary phase in the definition of U p z by making the choice 



(50) 

(51) 
(52) 
(53) 
(54) 

(55) 



It is readily confirmed that Uk and Up are order 2. However Ua, unlike A, is only 
order 12. This is because [25] the mapping G h-> Ug takes A 12 = ( § § ) to a multiple 
of the identity. 

We seek the joint eigenvectors of the commuting unitaries U P and U\. We have 



TI 2 



U A U A 



Moo 








°\ 











A 12 





A 23 








V o 





A 3 i 


o / 



where 



A rs — A r yl s 



(56) 



(57) 



The 2x2 block Aqq is easily diagonalized. To diagonalize the other, 6x6 diagonal 



block of U\ observe that 



A12A31A23 — A23A12A31 — A31A23A12 — I 



Consequently 




A u 
A 23 
A31 

with £ 7^ if and only if £ is a cube root of unity and 

Ul = CA12A31U2 U3 = t; 2 A 3 iU2 

If we impose the further requirement that 

Ul 
U2 



(58) 



(59) 



(60) 



(61) 
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be an eigenvector of Up we must have that U2 is a multiple of ( J ) or ( \ ) . Wc 



conclude that a complete set of joint eigenvectors is 

|0,l) = i=(0, 0, ry 33 , ^\ W 5 , 0, -1, l)* (62) 

|0,-l) = i(0, 0, V 27 , i] 3 , 0, V2v 33 , V 6 , V e Y (63) 
V6 

|1,1> = | (/7 23 ^/6(3 + V3) 5 ^^6(3 -V3), 0, 0, 0, 0, 0, o)' (64) 

|2,l} = -±=(0, 0, 77 17 , r? 17 , vV\ 0, -1, l)* (65) 

|2,-l) = -^(0, 0, rj 11 , V 35 , 0, V2 V , rf, rfY (66) 

|4,l) = -^(0, 0, 77, V , W 7 , 0, -1, 1)* (67) 

|4,-1) = -L(0 5 0, r, 43 , ^ 0, V2ry 17 , r? 6 , r? 6 )' (68) 

|5,1) = \ (77^6(3- V3), r/ 3 y6(3 + V3), 0, 0, 0, 0, 0, 0) (69) 
where 77 = e5i , 

C/i|r,s)=7 ? 8 '-|r, S ) C/p|r, S ) = S |r,s) (70) 
for all r, s, and where the phases have been chosen so that 

U K \r,s) = \r,s) (71) 



for all r, s. Note that the fact that Uk is an anti-unitary means that Eq. (|7Tj) fixes 
the phase of |r, s) up to a sign (this is an important point: if one chooses the phases 
differently the fiduciality conditions are much harder to solve) . The action of the 



anti-unitary U A is: 

U A \0,1) = 77 42 |0, 1> U A \0,-1) = 77 30 |0, — 1) (72) 

U A \1,1) =?7 14 |5,1) (73) 

U A \2,1) ^ 10 |4,1) C/a|2,-1) =?7 46 |4,-1) (74) 

C^|4,l> ^ 26 |2,1) C/ A |4,-l) =77 14 |2,-1) (75) 

U A \5,1) =77 22 |1,1) (76) 



The eigenspaces of U F — U\ are as in Table [T] As we remarked above N = 8 is 



eigenspace 


eigenvalue 


dimension 


basis 


So 


1 


2 


[0,1), [0,-1) 


Si 


r? 16 


3 


[5,1), [2,1), [2,-1) 


s 2 


r, 32 


3 


|1,1>, |4,1>, |4,-1) 



Table 1. 



one of the dimensions for which fiducials exist in all 3 eigenspaces of the Zauner 
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unitary. Since U A toggles the eigenspaces S\ and 52 we only need to calculate 
fiducials in So and Si . We thus need to solve the equations 

,2 1 



D r 



(77) 



'PIV/I — g 

for all p 7^ (mod 8) and \ip) of the form 

\tp) = cos0|5,l) +sin6»cos< ? !)e 4X2 |2, 1) + sin 9 sin <jye lX3 |2, -1) (78) 

with < 6, 4> < § and < %2, X3 < 2tt, or 

|V>) = cos0|O,l) + sin(9e JX |0,-l) (79) 

with < 9 < | and < x < 27r. The fact that \ip) is assumed to be an eigenvector 
of Uf z means that it is enough to require that it satisfy the 11 equations 

,2 1 



^\D V ^)\ = 
Da- P2 



\(i>\D A r p3 lP) 

\y,\D A r Pi i>) 



r = 0, 1 
r = 0,...,3 
r = 0,...,3 



with 



Pi = (0,4) P2 = (0,2) P3 = (0,l) P4 = (l,2) 

or, equivalently, the following linear combinations of these equations: 

\(iP\D P2 iP)\ 2 + \(iP\D Ap2 i;) ' 



r=0 



r=0 



r=0 r=0 



r=0 



r=0 



|(V>I^ P2 V) 

(V|^ A r p4 V>) 

(VI^ P4 V) 
(V|^ P4 V) 

(Vl^psV) 

(V|^ P4 V) 



(80) 
(81) 
(82) 
(83) 

(84) 
(El) 



= (E2) 

-§ (E3) 

= (E4) 

= (E5) 

= (E6) 

= (E7) 

= (E8) 

= (E9) 



r=0 
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^fcW|<#DA^>| 2 = (E10) 

r=0 
3 

]Tfc( 4 )|(V^p 4 V>}f = o (eh) 



where 



1 
1 

w 



£;( 2 > 



-1 
1 

v-v 



r=0 



-1 
V o / 



few 



1 



v-v 



(85) 



We begin with the fiducials in Si. Substituting Eq. ([78)1 into these equations we 
find that Eqs. (|Elj) - (|E3j) are each equivalent to the single condition 



cos 26>(1 + cos 26>) = 



(86) 



The solution 9 = ? is inconsistent with the remaining equations. On the other hand 
if one sets 9 = \ in Eqs. (|E4j) — ()E11[) it is straightforward (though perhaps a little 
tedious) to show that the resulting equations are equivalent to the four conditions 



sin 2 2</>cos2(x2 — Xa) = — 2-\/3cos20 

1 



(87) 



cos </>cos2x2 = 



sin 6 cos 2%3 = - 
cos 2 20-4cos20- 1 = 
whose solution is 



_ cos z 26 + 2 cos 26 - 1 

1/9 S 

cos 2 26- 2cos2<^)- l) 



4^/2 



COS = 



'3-V5 



,1X2 



A' 2 



8 + V6 - ViO + isi V 8 - \/6 



30 J 



5 ( A3 — p 12 

" 4 



V6 - VZO + isi 



V6- 



(89) 
(90) 

(91) 
(92) 
(93) 



where si, S2, S3 = ±1 can be chosen independently. 

Turning to the fiducials in Sq we find on substituting Eq. (|T9[) into Eqs. (|BT 
(jElip that the equations reduce to the two conditions 



sin 2 29 cos 26 = 
1 + cos 26» - cos 2 26> = 



whose solution is 



= e 



(94) 
(95) 

(96) 



where r — 0, 



.3. 
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To conclude: there are 8 fiducials in the 1S1 eigenspace given by 



r/ 37 \/3 



V3 





= 

2^/3 



\ o 

with si, S2, S3 = ±1 and 




isix 






v n 


-1 

V i / 




l(x-T5) 



(97) 



4 



VG- V36+i\/8- Ve + Vzd) (98) 



another 8 fiducials in the S2 eigenspace obtained by acting on these with the anti- 
unitary Ua, and 4 fiducials in the Sq eigenspace given by 




/ \ 


?7 33 

yy 33 

V2rj 15 



-1 

V i / 




/ \ 

o 

?y 33 
rf 


vV 9 



n 



12 



(99) 



V v 12 J 



where r = 0, . . . , 3. The fiducials in £>i, 52 are on Scott-Grassl [36] orbit 8a; the 
fiducials in Sq are on Scott-Grassl orbit 8b. 

Comparing these fc-nomial basis expressions with the standard basis expressions 
in Ref. [35] we see that the degree of simplification achieved is very considerable. 



5. Dimension 12 fiducial 

Before we present our new solution for dimension 28, we revisit dimenion N = 
12 = 2 2 x 3 to further demonstrate the advantage of the sparse representation of 
the Clifford group. 

Choosing a basis such that both X 6 and Z 6 are diagonal, the representation of 
the Clifford group will be 3-nomial. In this basis, the Weyl-Heisenberg group is 
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generated by 



X 






w 24 






































, ,23 

oj 24 






































L0 24 






































, ,15 
W 24 






































1 1 17 
ui 24 




















1 1 7 






























































































, ,11 

L0 24 






































/ , 17 

W 24 






































W 24 






































U>24 




















oj 24 


















(100) 



and 



/ 








oj 2 





V o 











^24 

















U>24 







U), 














1 , 17 
w 24 



^24 















, 17 

^24 

















,7 

















UJ 24 

















,23 
^24 











o \ 









^24 







o / 



(101) 



where u>24 denotes a primitive 24th root of unity. 

Like in the case of dimension N — 8 discussed in detail in the previous section, 
we consider a maximal Abelian subgroup of the normalizer Mz of the Zauner sym- 
metry Fz in the Clifford group when choosing the basis of the eigenspace of Fz of 
dimension 5. We find the non-normalized orthogonal basis {vq, v\, v%, 1)3, V4,} with 
the vectors 



v = (0,0,0,0,0,0,0,1,0,1,0,1) 
vi = (1,1,1,1,1,1,1,0,1,0,1,0) 
v 2 = (1, 1,ojI,ojI,uj 3 ,uj 3 ,uj 3 ,Q, 1,0,^3,0) 

v 3 = (0,0,0,0,0,0,0,1,0,^,0,^3) 

1> 4 = (1,^3,^3, 1,^3, U3 3 , 1, 0,^3,0,^3,0) 



(102) 
(103) 
(104) 
(105) 
(106) 



In the approach for computing a solution with Zauner's symmetry reported in 
[31j . computing a single modular Grobner basis took about 40 hours using more 
than 17 GB of memory. Using the Abnormal representation of the Clifford group 
and an adapted basis of the eigenspace of Zauner's matrix, the corresponding step 
uses less than 100 MB of memory and takes less than 40 seconds on the same 
hardware-reduced to less than 25 seconds on current computers. Moreover, instead 
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of more than 300 different primes, now only 21 suffice. Overall, we get a speed-up 
by more than 3 orders of magnitude. 
A non-normalized fiducial vector 

n 

|^o>=$>^> (107) 

z=0 

with respect to the representation of the Weyl-Heisenberg group generated by the 
matrices X and Z given in Eqs. (|100j) and (|101j) . respectively, can be expressed in 
the number field 

K = Q(V2,V3, VT3,t 1 ,s 1 ,s 2 ,V ZI T) (108) 

of degree 192, where 

si = \J{Vl3- l)/2, s 2 = ^(3^ + 9)72, and t\ = \2t x + 10. (109) 
In this field, a primitive 24th root of unity is given by 

uj 2A = il ((i - VS) + (V3 + i)V=T) ■ (no) 

The coefficients of \ijjo) in Eq. f| 107[) are as follows: 

x = ((-30v / 13 - 312)si + (24-v/l3tf - 102-\/l3ti - 309%/l3 - 92t? - 158ti - 5)) V / ^T 

+ (26\/39 - 364^)51 + 28^4? - 58^*1 - 147^ + 96^*? - 43V3ti - 443^ 

xi = (24v / 13<? - 102v / 13ii - 540v / 13 - 92t 2 L - 158ti - 980) V / ^T 

+ 28-v/39ti - 58^*1 - 264V39 + 96\/3ti - A2Vlt 1 - 1184\/3 

x 2 = ((24 v / 13 - 702)si - 54v / 13t? + 138Vl3ti + 375v / 13 - 98*? + 142ti + 667) V /T T 

+ (28v / 39 - 26^)^1 - 2v / 39tf - 22v / 39ii - 81%/39 - 94^*? - 58\/3ii + 219^ 

x 3 = ((-30v / 13 - 312)«i - 54^*1 + 138^*1 + 315^ - 98*? + 142ti + 43) v^T 

+ (26^ - 364^)31 - 2-v/39ti - 22v / 39ti + 93^ - 94 v / 3t? - 58\/3*i + 1077^ 

x' 4 = (30v / 13t? - 36Vl3ti - 588v / 13 + 1904? + 16ti - 3236) 

- 26-v/39t? + 80^*1 + 168^ - 2V~3t 2 + 10Ch/3ti - 400^ 

x 5 = ((24^- 702)81 + (30v / 13<I - 36-v/l3ti - 297v / 13 + 190i? + 16ti - 1637)) v^T 

+ (28V39 - 26^)^1 - 26V39t'( + SOv^ti + lllv / 39 - 2V~3t'( + 100"v/3ti - 517^3 

X6 = ((-30-/13- 312)si + (3Q-v/l3tf - 36^*1 - 357^+1904? + 16*i - 2261))^" 

+ (26^ - 364^)81 - 26^^ + 80^*1 + 285^ - 2v / 3t? + 10uV3ti + 341^ 

x 7 = 488v / 13s 2 

x$ = ((24Vl3 - 702)si + (24v / 13tf - 102\/l3ti - 249^ - 92t? - 158ti + 619))v / ^T 

+ (28v / 39 - 26^)^1 + 28v / 39t? - 58-^*1 - 321v / 39 + 96%/3t? - 43y/3ti - 1301-\/3 

Xg = ((85^ + 91^)8152 - 122v / 39s 2 )v /3 T 
+ (23\/l3- 871)sis 2 + 122v / 13s2 

xio = (-54-v/l3t? + 138-^*1 + 84VT3 - 98t? + 142ti - 932) v /:r T 

- 2v / 39i? - 22^*1 - 24^- 94\/3ii - 58\/3ii + 336^3 
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xii = ((3lV39 + 481\/3)sis 2 + 122V39s2)^ 
+ (139-/13 - 299)sis 2 + 122VT3S2 



(111) 



This expression for the fiducial vector is much more compact than the page filling 
ones given in Ref. .31] and the supplementary material of Ref. |36j . 

6. Dimension 28 fiducial 

In order to compute a fiducial vector for dimension N = 28 = 2 2 x 7 with the 
help of the computer algebra system Magma [70], we first compute a basis such 



that both X 14 and Z 14 are diagonal, 
given by the matrix 



It turns out that such a change of basis is 



T 



-(' 



/o 








1 









o\ 







1 








0/ 



(112) 



which is a tensor product of a 2 x 2 Fourier transformation and a permutation 
re-ordering the basis elements with odd/even index. The numerical solution in the 
orbit 28c in |36] indicates that there is a solution which possesses an order-two 
anti-unitary symmetry 

'11 6 



F K = 



50 17 



(113) 



in addition to the symmetry F z of order three conjectured by Zauner. In order 
to compute a basis of the corresponding eigenspace, we represent the anti-unitary 
transformation Uf c f z as a real orthogonal matrix Of c f z € 0(2N) using the isomor- 
phism C — > M 2 . We obtain ten linearly independent eigenvectors V t 6 M 56 which 
are then mapped back to ten vectors G C 28 . Then we write the fiducial vector 
as 



(114) 



i=0 



with real variables j/j, i = 0, . . . , 9. The defining equations for the SIC are given as 

(^ol^o) = 1 and (115) 



|(Vo|C P |^o)| 2 = 4> forp^O 



(116) 



Pirw, 2 g: 

We solved these equations via computing a Grobner basis for the corresponding 
polynomial ideal, using similar techniques to those described in |31) . First, the 
equations are reformulated as equations over the rationals, adding an auxiliary 
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variable. This allows us to compute a modular GB, i.e., performing all computations 
modulo a large prime. 

Using Magma V2.18 on a Linux PC with 3 GHz clock speed and 32 GB RAM, 
the first step took about 60 hours and used up to 45 GB of memory. One of the 
polynomials in the modular Grobner basis could be lifted to the rationals. Adding 
this polynomial, the computation for three more primes took only approximately 
six hours. Having computed four modular Grobner bases, six polynomials could be 
lifted to the rationals. Adding those polynomials reduced the time to compute a 
modular Grobner basis to only five minutes. Then, with five different primes, we 
lifted 20 polynomials to the rationals. The following Grobner basis computations 
took less than one second each. In total we used 18 different primes with about 23 
bits each to obtain a Grobner basis over the rationals. As a next step, the Grobner 
basis is converted to a Grobner basis with respect to lexicographic order. This 
yields a polynomial system of equations in triangular form. The coefficients of the 
polynomials have numerators and denominators with more than 50 digits. It turns 
out that the ideal is zero-dimensional, i.e., there are finitely many solutions. Note 
that we are only interested in solutions where all variables y\ assume real values 

Inserting one of these solution into Eq. (|114j) . we find that a non- normalized 
fiducial vector 

27 

|^o) = X>^> 

i=0 

in the adapted bases given by the matrix T in Eq. (|112j) can be expressed in the 
number field 

K = Q(v^,V7,V29 ) ri,r 2 ,«i ) fl 2 ,>/ = l) (118) 

of degree 576, where 

n = yV29 + 5, r 2 = \JV29 + l, sf = 21si-7, and s% = 14s 2 -14. (119) 

Note that even though we do not explicitly use a 56th root of unity, it can be 
expressed as 

^56 = ((«? - 7ai - 35) (1 + V=l)Vu - (7s\ - 7si - 77) (1 - V^1)V2) . (120) 

The Galois group of K is isomorphic to Cq x {{Cq x Ci x C 2 x C2) x Cg), and K is 
an Abelian extension of Q(y29). 

The coefficients of |-0 O ) m Eq. (|117[) are as follows: 

x = (((2a 2 «i + (4s 2 - 21)si - 7s 2 + 42)^29 

+ ((-12s 2 - 18s 2 + 112)s? + 203si + (210s^ + 315s 2 - 2366)))v / 2r 2 
+ ((l/7(4s§ + 28s 2 )si + (8s§ + 20s 2 - 84)si - 26s| - 110s 2 + 

+ (l/7(60s| - 84s 2 - 560)sf + (-24s| - 36s 2 + 224)si - 102s| + 282s 2 + 952))v / 7)ri v^T 

+ (((9*i + 9s 2 )v / 29 - 87si + 18^2 + 27s 2 - 168)v / 14r 2 

+ ((-84si + (18s| + 126s 2 - 336))^ + (270si - 378s 2 - 2520) ))n 

xi = (((-s 2 «i + (-2s 2 + 21)si + (35s 2 + 84))v / 29 

+ ((-6« 2 - 9s 2 + 56)s? - 203si + 42s| + 63s 2 - 1204))V2r 2 

+ ((1/7(10«| + 28s 2 - 168)sf + (8s| + 26s 2 - 84)si - 14s| - 68s 2 + 336)\/29 

+ (l/7(-18s| - 56s 2 + 168)«i - 58s 2 si - 90s! + 184s 2 + 840))v / 7>i ^1 
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+ ((((s 2 - 2)s\ + (2s 2 - l)si - 17s 2 + 28)^ 

+ ((-2s! - 3s 2 + 38)s? + (8s! + 12s 2 - 65)si + (10s! + 15s 2 - 364)))v / 14r 2 

+ (((- 2s 2 - 8s 2 )sf + (8s| + 14s 2 - 84)si + (46sl + 112s 2 - 336))v / 29 

+ ((-6s! + 20s 2 + 56)s? + (-48s! - 14s 2 + 448)si - 150s| - 196s 2 + 1400)))n 

x 2 = ((2s 2 s? - 8s 2 si - 28s 2 )v / 58r 2 

+ ((l/7(8aa + 14s 2 - 112)s? + (4s| + 16s 2 - 56)si - 28s! - 40s 2 + 392)v / 29 

+ (l/7(-48sa - 14s 2 + 448)s? + (12s! - 40s 2 - 112)si + (204s! + 16s 2 - 1904))) \/7)n V^l 

+ (((-2s? - 431+70)^ 

+ ((4s| + 6s 2 - 18)s? + (8s| + 12s 2 - 36)si - 32s! - 48s 2 - 378))v / 14r 2 

+ (((- 4s 2 - 10 *2 + 56)sf + (4s| + 28s 2 - 56)si + (116s! + 308s 2 - 1120))v / 29 

+ ((12s! + 18s 2 - 112)si + (60sl - 84s 2 - 560)si - 276s! - 588s 2 + 2576)))n 

x 3 = ((l/7(-16s| - 28s 2 + 168)s? + (4s| + 16s 2 )si + (68s! + 128s 2 - 672))v / 29 

+ (l/7(96s| + 28s 2 - 896)s? + (12s! - 40s 2 - 112)si - 372s! - 152s 2 + 3472))A/7riA/^T 
+ (((-§ s 2 - 20S2 + 56)s? + (-52s! - 112s 2 + 448)si + (172s! + 448s 2 - 1120))v / 29 
+ ((24s! + 36s 2 - 224)s? + (228s! + 168s 2 - 2128)si - 444s! - 840s 2 + 4144))n 

x 4 = ((-2s 2 s? + 8s 2 S! + 28s 2 )v / 58r 2 

+ ((V 7 ( 8s 2 + 14s 2 - 112)s? + (4s! + 16s 2 - 56)si - 28s! - 40s 2 + 392)v / 29 

+ (l/7(-48s| - 14s 2 + 448)s? + (12s! - 40s 2 - 112)si + (204s! + 16s 2 - 1904)))v / 7)n\/^l 

+ (((2sf +4si - 70)^ 

+ ((-4sl - 6s 2 + 18)s? + (-8s! - 12s 2 + 36)si + (32s! + 48s 2 + 378)))v / 14r 2 
+ (((- 4s ! - 10 *2 + 56)s? + (4s! + 28s 2 ~ 56)si + (116s! + 308s 2 - 1120))\/29 
+ ((12s| + 18s 2 - 112)s? + (60s! - 84s 2 - 560)si - 276s! - 588s 2 + 2576)))n 

x 5 = (((s 2 sf + (2s 2 - 21)si - 35s 2 - 84)v / 29 

+ ((6s! + 9s 2 - 56)s? + 203si - 42s! - 63s 2 + 1204))v / 2r 2 
+ ((l/7(10s! + 28s 2 - 168)s? + (8s! + 26s 2 - 84)si - 14s! ~ 68s 2 + 336)729 
+ (l/7(-18s! - 56s 2 + 168)ai - 58s 2 si - 90s! + 184s 2 + 840))\/7>i V^l 
+ ((((-s 2 + 2)s2 + (-2s 2 + l)si + (17s 2 - 28))%/29 

+ ((2s! + 3s2 - 38)s? + (-8s! - 12s 2 + 65)si - 10s! - 15s 2 + 364))Vl4r 2 

+ (((-2s! ~ 8s 2 )sf + (8s! + 14s 2 - 84)si + (46s! + 112s 2 - 336))v / 29 

+ ((-6s! + 20s 2 + 56)s? + (-48s! - 14s 2 + 448)si - 150s! - 196s 2 + 1400)))n 

x 6 = (((-2s 2 s? + (-4s 2 + 21)si + (7s 2 - 42))^29 

+ ((12s! + 18s 2 - H2)sf - 203si - 210s! - 315s 2 + 2366))v / 2r 2 

+ ((!/7(4s! + 28s 2 )sf + (8s! + 20s 2 - 84)si - 26s! - 110s 2 + 168)^29 

+ (l/7(60s! - 84s 2 - 560)s? + (-24s! - 36s 2 + 224)si - 102s! + 282s 2 + 952))V7y ri V^l 

+ (((-9si - 932)^ + (87si - 18s! - 27s 2 + 168))v / 14r 2 

+ ((-84si + (18s! + 126s 2 - 336))^ + (270s! - 378s 2 - 2520)))n 



(((l/7(10s| + 28s 2 - S4)sl + (2s| + 2s 2 )si - 38s! - H0s 2 + 336)^ 
+ (l/7(-18s| - 56s 2 + 168)s? + (-18s| + 2s 2 + 168)si + (54s! + 226s 2 - 504)))v / 7 
((-2s! - § s 2 + 28)sf + (-22s! - 70s 2 + 224)sj + (34s! + 154s 2 - 560))\/29 

-6s| + 20s 2 + 56)s? + (6s| + 154s 2 - 56) si + (174s| - 406s 2 - 1624))))a/2>ia/ = 
((l/7(10s| + 28s 2 - 84)s? + (2s| + 2s 2 )si - 38s! - 110s 2 + 336)\/29 
(l/7(-18s| - 56s 2 + 168)sf + (-18s| + 2s 2 + 168)si + (54s! + 226s 2 - 504)))v / 7 
((2s! + 8s 2 - 28)sf + (22s! + 70s 2 - 224)si - 34s! - 154s 2 + 560)V29 
(6s| - 20s 2 - 56)s? + (-6s! - 154s 2 + 56)si - 174s| + 406s 2 + 1624)))v / 2r 1 



((( 



((( 



s 2 - 2) si + (2s 2 - 4)si - 8s 2 + 70)^ 

-2s! - 3s 2 + 38)s^ + (-4s! - 6s 2 + 76)si + (16s! + 24s 2 - 826)))\/7 
s 2 Sj - 4s 2 si - 14s 2 )v / 29 

6s| + 9s 2 - 56)s'f + (-24s! - 36s 2 + 224)si - 84s! - 126s 2 + 784)))r 2 
(l/7(-5s! - 14s 2 + 28)s( + (2s| + 2s 2 - 28)si + (22s! + 58s 2 - 140))a/29 
l/7(9s| + 28s 2 - 84)s^ + (-18s! + 2s 2 + 168)si - 54s! - 110s 2 + 504))^ 
(-s| - 4s 2 )s^ + (-8s! - 14s 2 + 84)si + (20s! + 98s 2 - 168))a/29 
-3s! + 10s 2 + 28)s? + (48s! + 14s 2 - 448)si + (132s! - 266s 2 - 1232))))v / 2)n 1 
(((-s 2 + 2)sl + (-2s 2 + 4)si + (8s 2 - 70))\/29 
2s! + 3s 2 - 38)s? + (4s! + 6s 2 - 76) si - 16s! - 24s 2 + 826))^ 
S2sj - 4s 2 si - 14s 2 )v / 29 

6s! + 9s 2 - 56)sf + (-24s! ~ 36s 2 + 224)si - 84s! - 126s 2 + 784)))r 2 
(l/7(-5s| - 14s 2 + 28)s^ + (2s! + 2s 2 - 28)si + (22s! + 58s 2 - 140))a/29 
l/7(9s| + 28s 2 - 84)s^ + (-18s! + 2s 2 + 168)si - 54s! - 110s 2 + 504))^ 
(s! + 4s 2 )s? + (8s! + 14s 2 - 84)si - 20s! - 98s 2 + 168)a/29 
3s! - 10s 2 - 28)s^ + (-48s! - 14s 2 + 448)si - 132s! + 266s 2 + 1232))) v^n 



-s 2 - 2)sl + (s 2 - l)si + (lls 2 + 28))\/29 

-2s! - 3s 2 + 38)s^ + (-10s! - 15s 2 + 103)si + (46s! + 69s 2 - 700)))a/7 
S2s\ + (-82 - 21)si + (7s 2 - 84))\/29 

-6s! - 9s 2 + 56)s^ + (-6s! - 9s 2 + 259)si + (126s! + 189s 2 - 364))))r 2 
(l/7(-s| - 7s 2 - 28)s? + (s| - 2s 2 - 14)si - 13s! - 19s 2 + 224)\/29 
l/7(-15s| + 21s 2 + 140)sf + (-21s! + 12 S2 + 196)si + (93s! + 9s 2 - 868)))a/7 
(si + s 2 )s! + (5s| + 14s 2 - 42)si + (13s! + 49s 2 - 168))V29 
-9s! + s 2 + 84)s? + (-9s! - 28s 2 + 84)si + (27s! - 119s 2 - 252))))v / 2)ri 
(((s 2 + 2)s\ + (-82 + l)si - lls 2 - 28)\/29 

2s! + 3s 2 - 38)s? + (10s! + 15s 2 - 103)si - 46s! - 69s 2 + 700)) 
S2s'i + (-82 - 21)si + (7s 2 - 84))\/29 

-6s! - 9s 2 + 56)s^ + (-6s! - 9s 2 + 259)si + (126s! + 189s 2 - 364))))r 2 
(l/7(-s! - 7s 2 - 28)sf + (s| - 2s 2 - 14)si - 13s! - 19s 2 + 224)\/29 
l/7(-15s| + 21s 2 + 140)s? + (-21s! + 12s 2 + 196)si + (93s! + 9s 2 - 868)))^7 
((-si - s 2 )s? + (-5s! - 14s 2 + 42)si - 13s! - 49s 2 + 168)\/29 
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- ((9s! - s 2 - 84)s? + (9s? 2 + 28s 2 - 84) si - 27s! + 119s 2 + 252)))V2)n 



xio = (((9sia/29- 87si - 36s! - 54s 2 + 336)\/7 + ((4s 2 sf + (2s 2 + 21)si - 56s 2 - 42)729 
(12s! + 18s 2 - 315)si + (84sl + 126s 2 - 378))))r 2 

((l/7(8sl + 14s 2 - 56)4 + (-2*1 - 8«2 + 14)si - 16s! - 19s 2 + 112)v / 29 
l/7(-48sa - 14s 2 + 448)s? + (-6s! + 20s 2 + 56) si + (132s! - 5s 2 - 1232)))\/7 
42si - 36s! - 63s 2 + 420)\/29 + (216s| + 63s 2 - 2016)))v / 2)ri v^T 



xii = ((( 



+ 



X12 



((( 



(-9sia/29 + (87si + (36s! + 54s 2 - 336)))\/7 + ((4s 2 s? + (2s 2 + 21)si - 56s 2 - 42)\/29 
12s| + 18s 2 - 315)si + (84s! + 126s 2 - 378))))r 2 

(l/7(8sl + 14s 2 - 56)s? + (-2s! - 8s 2 + 14)si - 16s! - 19s 2 + 112)v / 29 
l/7(-48s| - 14s 2 + 448)s? + (-6s| + 20s 2 + 56) si + (132s! - 5s 2 - 1232)))v / 7 
42si + (36s! + 63s 2 - 420))\/29 + (-216s! - 63s 2 + 2016)))\/2>i 



9si\/29 + (87si + (36a! + 54s 2 - 336)))a/7 + ((-4s 2 s? + (-2a 2 - 21)si + (56s 2 +42))V29 
-12s! - !8s2 + 315)si - 84s! - 126s 2 + 378)))r 2 

(l/7(8s! + 14s 2 - 56)s? + (-2s| - § s 2 + 14)si - 16s! - !9s 2 + 112)%/29 
l/7(-48s! - 14s 2 + 448)s? + (-6s! + 20s 2 + 56) si + (132s! - 5s 2 - 1232)))\/7 
42si - 36s! - 63s 2 +420) a/29 + (216s! + 63s 2 - 2016)))v / 2>iv /r T 
(9si\/29 - 87si - 36s! - 54s 2 + 336)\/7 + ((-4s 2 s? + (-2s 2 - 21)si + (56s 2 + 42))\/29 
-12s! - 18s 2 + 315)si - 84s! - 126s 2 + 378)))r 2 

(l/7(8s! + 14s 2 - 56)s? + (-2s! - 8*2 + 14)si - 16s! - 19s 2 + 112)a/29 
l/7(-48s! - 14s 2 + 448)s? + (-6s! + 20s 2 + 56) si + (132s! - 5s 2 - 1232)))v / 7 
42si + (36s! + 63s 2 - 420))\/29 + (-216s! - 63s 2 + 2016)))\/2>i 

s 2 + 2)s^ + (-S2 + l)si - lls 2 - 28)\/29 

2s! + 3s 2 - 38)s? + (10s! + 15s 2 - 103)si - 46s! - 69s 2 + 700))\/7 
-s 2 s? + (s 2 + 21)si - 7s 2 + 84)a/29 

6s! + 9s 2 - 56)s? + (6s! + 9s 2 - 259)si - 126s! - 189s 2 + 364)))r 2 
(l/7(-s! - 7s 2 - 28)s? + (s! - 2s 2 - 14)si - 13s! - 19s 2 + 224)\/29 
l/7(-15s! + 21s 2 + 140)sj + (-21s! + 12s 2 + 196)si + (93s! + 9s 2 - 868)))a/7 
(s! + s 2 )sl + (5s! + 14s 2 - 42)si + (13s! + 49s 2 - 168)) V29 
-9s! + s 2 + 84)s? + (-9s! - 28s 2 + 84)si + (27s! - 119s 2 - 252)))) v / 2)ri\/^I 
(((-s 2 - 2)s\ + (s 2 - l)si + (lls 2 + 28))a/29 

-2s! - 3s 2 + 38)s^ + (-10s! - 15s 2 + 103)si + (46s! + 69s 2 - 700)))a/7 
-s 2 sf + (s 2 + 21)si - 7s 2 + 84)\/29 

6s! + 9s 2 - 56)s? + (6s! + 9s 2 - 259)s x - 126s! - 189s 2 + 364)))r 2 
(l/7(-s! - 7s 2 - 28)s? + (s! - 2s 2 - 14)si - 13s! - 19s 2 + 224)\/29 
l/7(-15s! + 21s 2 + 140)s? + (-21s! + 12s 2 + 196)si + (93s! + 9s 2 - 868)))\/7 
{-si - s 2 )s'f + (~5s! - 14s 2 + 42)si - 13s! - 49s 2 + 168)\/29 
9s! - si- 84)s? + (9s! + 28s 2 - 84) si - 27s! + 119s 2 + 252)))v / 2>i 
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sis = (((((-S2 + 2)s? + (-2S2 + 4)si + (8s 2 - 70))\/29 

+ ((2s! + 3s 2 - 38)s? + (4s! + 6s 2 - 76) si - 16s! - 24s 2 + 826)) + ((-s 2 s? + 4s 2 si + 14^2)^ 

+ ((-6«2 - 9s2 + 56)s? + (24s! + 36s 2 - 224)si + (84s! + 126s 2 - 784))))r 2 

+ (((l/7(-5s| - 14s 2 + 28)s? + (2s| + 2s 2 - 28)si + (22s! + 58s 2 - 140))V29 

+ (V 7 ( 9s 2 + 28s 2 - 84)s? + (-18s! + 2s 2 + 168)si - 54s! - 110s 2 + 504))^ 

+ (((-si - 4s 2 )s? + (-8s! - 14s 2 + 84)si + (20s| + 98s 2 - 168))\/29 

+ ((-3s! + 10s 2 + 28)s? + (48s! + 14s 2 - 448)si + (132s! ~ 2 66s 2 - 1232))))v / 2>i^I 

+ (((((s 2 - 2)s? + (2s 2 - 4)si - 8s 2 + 70)^ + ((-2s! - 3s 2 + 38)s? 

+ (-4s! - 6s 2 + 76)si + (16s! + 24s 2 - 826))) ^7 + ((-s 2 s? + 4s 2 si + 1452)^ 

+ ((-6«! ~ 9s 2 + 56)s? + (24s! + 36s 2 - 224)sj + (84s! + 126s 2 - 784))))r 2 

+ (((l/7(-5s| - 14s 2 + 28)s? + (2s! + 2s 2 - 28)si + (22s| + 58s 2 - 140))729 

+ (l/7(9s! + 28s 2 - 84)s? + (-18s! + 2s 2 + 168)si - 54s! - 110s 2 + 504))^ 

+ (((s! + 4s 2 )s? + (8s! + 14s 2 - 84)si - 20s! - 98s 2 + 168)v / 29 

+ ((3s! - 10s 2 - 28)s? + (-48s! - 14s 2 + 448)si - 132s! + 266s 2 + 1232)))v / 2)ri 

x 14 = ((-lSsjv 7 ^ + (42si - 84)V29)y/2r 2 + (((16s? - 34si - 248)^29 - 174si - 696)^ 
+ ((-8s? - 322si - 476)^ - 232s? - 2030si + 812)))-/=T 

+ (-18s 1 v / 203 + (-42si + 84)v / 29)v / 2r 2 + ((-16s? + 34s x + 248)v / 29 + (174si +696))v / 7 
+ (-8s? - 322si - 476)v / 29 - 232s? - 2030si + 812 

xi5 = (((-4s? - 2si + 56)v / 203 + (42si + 168)v / 29)v / 2r 2 

+ (((40s? + 26sj - 788)729 + (232s? - 58si - 4988))^ 
+ ((60s? + 42si - 252) v / 29 + (116s? + 406si + 812))))-/+L 
+ ((-4s? - 2si + 56)v / 203 + (-42si - 168)v / 29)v / 2r 2 
+ ((-40s? - 26si + 788)^ - 232s? + 58s x + 4988)\/7 
+ (60s? + 42si - 252) V29 + 116s? + 406si + 812 

x 16 = ((4s? + 8s i - 140)v / 406r 2 + (((-16s? - 68si + 452)^- 34851 + 1740)\/7 + ((60s? - 84si)\/29 
+ (116s? - 812si + 3248))) V~ + (4s? + 8si - 140) v / 406r 2 

+ ((16s? +68si - 452)^29 + (348si - 1740) )%/7 + (60s? - 8451)^+1165? - 812si + 3248 

x 17 = (((-12s? - 144si + 192)a/29 - 116s? - 928si + 2320) \H + ((68s? + 112si - 1792)^ 
+ (348s? - 9744))) V / ^T+ ((12s? + 144si - 192)^29 + (116s? + 928si - 2320))\/7 
+ (68s? + 112si - 1792) \/29 + 348s? - 9744 

x 18 = ((-4s? - 8si + 140)v / 406r 2 + (((-16s? - 68si + 452)v / 29- 348si + 1740)^ 
+ ((60s? - 84si)v / 29 + (116s? - 812si + 3248))))^T 

+ (-4s? - 8si + 140) + ((16s? + 68si - 452)v / 29 + (348si - 1740))^ 
+ (60s? - 84si)\/29 + 116s? - 812si +3248 
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xig = (((4s? + 2si - 56)v / 203 + (-42si - 168)v / 29)v / 2r 2 

+ (((40s? + 26si - 788) V29 + (232s? - 58si - 4988)) \fi + ((60s? + 42si - 252)v / 29 
+ (116s? +406si + 812))))-/^!+ ((4s? + 2si - 56)v / 203 + (42si + 168)a/29)a/2> 2 
+ ((-40s? - 26si + 788)^29 - 232s? + 58s x + 4988) W 
+ (60s? + 42si - 252) V29 + 116s? + 406si + 812 

x 20 = ((18siv / 203 + (-42si + 84)\/29)a/2V 2 + (((16s? - 34si - 248)^29- 174si - 696)^ 
+ ((-8s? - 322si - 476)^ - 232s? - 2030si + 812)))-/=4 

+ (ISsiv 7 ^ + (42si - 84)v / 29)v / 2r 2 + ((-16s? + 34si + 248)v / 29 + (174si + 696))^ 
+ (-8s? - 322si - 476)v / 29 - 232s? - 2030si + 812 

X21 = (((-4s| - 4s 2 + 56)s? + (-8s! + 28s 2 + H2)si + (104s! + 140s 2 - 1456))a/29 

+ ((36s| - 4s 2 - 336)s? + (216s! - 140s 2 - 2016)si - 792s! - 28s 2 + 7392))riv /r T 

+ ((l/7(-4s| - 28s 2 + 56)s? + (-8s| - 20s 2 + 112)si + (8s! + 92s 2 - 112))v / 29 

+ (l/7(-60s| + 84s 2 + 560)s? + (24s! + 36s 2 - 224)si + (264s! - 300s 2 - 2464)))%/7Vi 

X22 = ((((s 2 + 2)s? + (2s 2 + 4)si - 8s 2 - 70)^ 

+ ((2s| + 3s 2 - 38)s? + (4s! + 6s 2 - 76) si - 16s! - 24s 2 + 826))v / 14r 2 

+ (((2s! + 2s 2 )s? + (-20s! - 56s 2 + 168)si - 76s! - 112s 2 + 672)v / 29 

+ ((-18s! + 2s 2 + 168)s? + (36s! + 112s 2 - 336)si + (540s! + 56s 2 - 5040))))r n /-L 

+ (((-s 2 s? + 4s 2 si + 14s 2 )v / 29 

+ ((6s! + 9s 2 - 56)s? + (-24s! - 36s 2 + 224)si - 84s! - 126s 2 + 784))\/2V 2 

+ ((l/7(-2s! - 14s 2 + 56)s? + (8s! + 20s 2 - 56) si + (16s! + 76s 2 - 280))\/29 

+ (l/7(-30s| + 42s 2 + 280)s? + (-24s! - 36s 2 + 224)s x + (96s! - 204s 2 - 896)))%/7)ri 

X23 = (((-2s? + (-3s 2 - l)si + (6s 2 + 28))\/29 

+ ((4s! + 6s 2 - 18)s? + (2s! + 3s 2 - 9)si - 56s! - 84s 2 + 252))v / 14V 2 

+ (((- 4s 2 - !0s 2 + 56)s? + (-2s! - 14 *2 + 28)si + (20s! + 14s 2 - 280))v / 29 

+ ((12s! + !8 S 2 - 112)s? + (-30s! + 42s 2 + 280)si - 204s! + 42s 2 + 1904)))n v^T 

+ (((-2s 2 s? + (-S2 - 21)si + (28s 2 - 84))v / 29 

+ ((6s! + 9s 2 + 147)si - 84s! - 126s 2 + 1596))v / 2r 2 

+ ((l/7(-8s| - 14s 2 )s? + (-10s! - 22s 2 + 84) si + (40s! + 106s 2 - 168))729 

+ (l/7(48s! + 14s 2 - 448)s? + (42s! + 34s 2 - 392)si - 96s! - 202s 2 + 896))v / 7)ri 

X24 = (((-9si +952)^+ (87si - 18s! - 27s 2 + 168))v / 14r 2 

+ ((-84si - 90s! - 252s 2 + 672) V29 + (162s! + 504s 2 - 1512)))n V^T 
+ (((2s 2 s? + (-2s 2 - 21)si - 49s 2 + 42)^ 

+ ((12s! + !8 S 2 - 112)s? + (12s! + i8s2 + 91)si - 126s! - 189s 2 + 770))i/2r 2 

+ ((l/7(20s| + 56s 2 - 224)s? + (4s! + 4s 2 - 28)si - 58s! - 148s 2 + 616)\/29 

+ (l/7(-36s| - 112S2 + 336)s? + (-36s! + 4s 2 + 336)si + (162s! + 272s 2 - 1512)))v / 7)r 1 
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z 25 = (((-9si +952)^+ (87si - 18s! - 27s 2 + 168))v / 14r 2 

+ ((84si + (90s! + 252s 2 - 672)) v / 29 + (-162s! - 504s 2 + 1512)))n v /:: T 
+ (((2s 2 Si + (-2s 2 - 21)si - 49s 2 + 42)^ 

+ ((12s! + 18s 2 - 112)si + (12s! + 18s 2 + 91)si - 126s! - 189s 2 + 770))v / 2r 2 

+ ((l/7(-20s| - 56s 2 + 224)s? + (-4s! - 4s 2 + 28)si + (58s! + 148s 2 - 616))\/29 

+ (l/7(36s! + 112s 2 - 336)s? + (36s! - 4s 2 - 336)si - 162s! - 272s 2 + 1512))v / 7)ri 

2:26 = (((-2s? + (-3s 2 - l)si + (6s 2 + 28))\/29 

+ ((4s! + 6s 2 - 18)s? + (2s! + 3s 2 - 9)si - 56s! - 84s 2 + 252))v / T4r 2 

+ ((( 4 *2 + !0S2 - 56)si + (2s| + 14s 2 - 28) si - 20s| - 14s 2 + 280)^29 

+ ((-12s! - !8 S 2 + 112)si + (30s! - 42s 2 - 280)si + (204s! - 42s 2 - 1904))))ni/-L 

+ (((-2s 2 s? + (-s 2 - 21)si + (28s 2 - 84))v / 29 

+ ((6s| + 9s 2 + 147)si - 84s! - 126s 2 + 1596))v / 2r 2 

+ ((l/7(8«2 + 14s 2 )sf + (10s! + 22s 2 - 84) si - 40s! - 106s 2 + 168)a/29 

+ (l/7(-48s| - 14s 2 + 448)s? + (-42s| - 34s 2 + 392)si + (96s! + 202s 2 - 896)))\/7>i 

X27 = ((((S2 + 2)s? + (2s 2 + 4)si - 8s 2 - 70)^ 

+ ((2s! + 3s 2 - 38)sf + (4s! + 6s 2 - 76)si - 16s! - 24s 2 + 826))VTir 2 

+ (((- 2s 2 - 2s 2 )sf + (20s! + 56s 2 - 168)si + (76s! + 112s 2 - 672))v / 29 

+ ((18s! - 2s 2 - 168)s? + (-36s! - 112s 2 + 336)si - 540s! - 56s 2 + 5040)))n v^T 

+ (((-s 2 s? + 4s 2 si + Msajv^ 

+ ((6s! + 9s 2 - 56)s? + (-24s! - 36s 2 + 224)si - 84s! - 126s 2 + 784))v / 2r 2 
+ ((l/7(2s! + 14s 2 - 56)s? + (-8s! - 20s 2 + 56) si - 16s! - 76s 2 + 280)v / 29 
+ (l/7(30s! - 42s 2 - 280)sf + (24s! + 30 S2 - 224)si - 96s! + 204s 2 + 896))v / 7)ri 

Finally, the (non-normalized) fiducial vector with respect to the standard basis can 
be obtained as = T\^!q). 

We have also tried to compute a fiducial vector for dimension iV = 18 = 3 2 x2. 
In that case, an even more sparse representation of the Clifford group with at most 
two entries per row and column can be obtained. The problem is that for N = 18 
we only have the Zauner symmetry (as indicated by the numerical solutions) , which 
means we have to search for a fiducial vector within a 7-dimensional complex vector 
space, resulting in 13 real parameters. By contrast for TV = 28 the additional 
symmetry reduces the number of real parameters to only 10. It turns out that the 
disadvantage, of having more parameters for N — 18, outweighs the advantage, 
of a more sparse group representation. Consequently we were unable to obtain a 
solution in the time available. 



In this paper we showed that the standard representation of the Clifford group 
has a system of imprimitivity consisting of fc-dimcnsional subspaces, where k is the 
square free part of the dimension. This means that one can choose a basis such that 
the representation matrices are all what we call fc-nomial, with exactly one non-zero 
k x k block in each row and each column of blocks. We then used this basis to 
perform a hand-calculation of the (already known) exact fiducials in dimension 8, 



7. Conclusion 
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and to obtain by machine-calculation a new exact fiducial in dimension 28. We 
also revisited the machine calculation of a fiducial in dimension 12 having only the 
Zauner symmetry, and found that the computation time was reduced by a factor 
of more than 10 3 , and that the solution obtained was much more compact. Our 
results suggest that the fc-nomial basis is better adapted to the SIC-problem than 
the standard basis. As we remarked in the Introduction, the fact that SICs exist 
at all (analytically or numerically) in every dimension < 67 in spite of the defining 
equations being greatly over-determined, and the fact that the Galois group of 
the known exact fiducials has a particularly simple form even among the class of 
solvable groups, both hint at the presence of some underlying symmetry or other 
algebraic feature of the equations which has hitherto escaped notice. Our hope is 
that the work reported here will help us to discover that hidden feature and so lead 
to a solution to the SIC-existence problem. 

The Clifford group has numerous applications. It appears to us that the sparsity 
of the representation matrices described here means that they are likely to be 
relevant to other problems, apart from the SIC-existence problem. 

As we mentioned in the Introduction there are several different, though closely 
related constructions which go by the name "Clifford group" . Besides the version 
of the group considered here there is the version defined to be the normalizer of the 
tensor product of several copies of the Weyl-Heisenberg group [TJ[2], anc ^ there is 
also what might be called the Galoisian version defined in prime power dimension 
using a finite field [3l|4] . It would be interesting to try to extend the analysis of this 
paper to this more general setting. 
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